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In both the Sun and the early universe, the 3 He(a,7) 7 Be reaction plays a key role. The rate of this 
reaction is the least certain nuclear input needed to calculate both the primordial 7 Li abundance 
in big bang nucleosynthesis (BBN) and the solar neutrino flux. Taking advantage of several recent 
highly precise experiments, we analyse modern 3 He(a,7) 7 Be data using a robust and minimally 
model dependent approach capable of handling discrepant data sets dominated by systematic rather 
O ' than statistical errors. We find 534(0) = 0.580 ± 0.043(0.054) keV b at the 68.3(95.4)% confidence 

O ' level. 
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I. INTRODUCTION 

^ ' As the nearest star, the Sun is the best studied. Models predict the neutrino fluxes produced by the radioactive 
D , decay of 7 Be and 8 B in the solar core Over the last decade, various neutrino observatories have measured the 
-H ' flux and flavour composition of the 8 B neutrinos from the Sun 0, Q , allowing constraints to be placed on neutrino 
mixing angles and mass-squared differences. New experiments sensitive to the 7 Be flux Q are ongoing. Both the 7 Be 
and 8 B neutrino flux predictions are nearly directly proportional to the astrophysical S factor for the 3 He(a,j) 7 Be 
reaction, S 34 , at a relative energy of ~ 20 keV; in fact, </>„( 7 Be) oc S , 34 (0) - 86 and <M 8 B) a ^(O) ' 81 UM- This 
sensitivity allowed S34 to be constrained using the measured solar neutrino fluxes and laboratory measurements of 
S17, the astrophysical S factor for the 7 Be(p, 7) 8 B reaction [?}■ The 3 He(a,7) 7 Be reaction is important not only in 
determining the solar neutrino fluxes, but also in other astrophysical environs. 

Primordial nucleosynthesis has been a robust prediction of hot big bang cosmology for over 40 years [1, [^, [l(| El, E3] • 
£s^j \ The theory explains the large universal abundance of 4 He as well as the origin of trace quantities of the light isotopes 
D, 3 He, and 7 Li. It is a theory with three free parameters, the cosmic baryon density, the neutron mean lifetime, and 
I the number of active light neutrino species. Via fits to the standard electroweak theory, measurements at the Large 
Electron Positron Collider have determined the number of active light neutrino species to be N v = 2.984 ± 0.008 [l3j, 
justifying the choice N v = 3. Many experiments have been performed to determine the mean lifetime of the neutron; 
the current Particle Data Group recommendation is r„ = 885.7 ± 0.8 sec [l3j |. The cosmic baryon density has been 
determined by analysing anisotropies in the cosmic microwave background radiation; the latest WMAP satellite 
results are Vl^h 2 — 0.02273 ± 0.00062 [3], where Sl^ is the universal mean density of baryons expressed in units 
of the critical density required to close the universe and the Hubble constant is 100ft- km sec -1 Mpc" 1 . Using the 
Friedmann- Robertson- Walker cosmological model, the standard model of particle physics, and nuclear reaction rates 
one can therefore predict the light element abundances very precisely and compare directly with primordial abundance 
estimates based on observations of the oldest systems. 

When observations of the Li abundances in the oldest stars in the Milky Way were compared with the predictions 
of BBN, a 2-3cr discrepancy was found [1 51 ] . Several possibilities exist for resolving this discrepancy, including particle 
physics beyond the standard model, improvements to the stellar models used to interpret astronomical observations, 
better approximations of the matter distribution in relativistic cosmology [iH . fl7| . modifications of gravitational 
theory (e.g., [Hj]), and improved nuclear reaction rates. In this paper, we examine the last possibility. Lithium is 
made as beryllium in the early universe via the 3 He(ev,7) 7 Be reaction. Of course the rates of reactions that destroy 
7 Be must be known in addition to those that create it. Recent studies of the 7 Be(c?, p)2a reaction cross section suggest 
that this reaction is not the source of the discrepancy [l^, [2fj| . The primordial 7 Li abundance prediction is nearly 
directly proportional to the 3 He(a,7) 7 Be cross section at a relative energy of ~ 300 keV; the primordial abundance 
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ratio 7 Li/H oc S®f 6 [l5| . On account of the obvious importance of this reaction for both solar neutrinos and BBN, it 
has been subjected to extensive theoretical as well as experimental study. 

many theoretical studies of the 3 He(a,7) 7 Be reaction have 



Since the first potential model calculations 
been performed [|| H US S S3, SI SH M 



21, U 

31] . Examining the potential models, a hard sphere potential yields 
5 34 (0) = 0.47 keV b and S"(0)/5(0) = -0.60 MeV" 1 (H, while a more physical potential yields 634(0) = 0.46 keV b 
and S'(0)/S(0) = —0.79 MeV -1 31]. Though the zero-energy S factor values are in good agreement, the shapes as 
measured by the logarithmic derivatives are quite different. Similarly, in the single channel resonating group method 
(RGM) calculations of Ref. [25] , 5(0) varies widely with the nucleon-nucleon potential. The most commonly used 
fitting function for experimental data is the microscopic cluster model calculation of Kajino [H SE S3]- Ref. (25j 
reports 634(0) = 0.50±0.03 keV b and S'(0)/S(0) = -0.548±0.033 MeV" 1 . The estimated theoretical uncertainty of 
±6% [25[ on both of these quantities must be considered when fitting this predicted shape to experimental data. The 
commonly cited potential model and microscopic cluster model calculations are shown in Fig ure [T] Another RGM 
calculation that includes the 6 Li+p configuration in addition to the 3 He+ct configuration [30] finds that the energy 
dependence of the astrophysical 6 factor is not uniquely determined and yields a range of 6'(0)/6(0) from —0.70 to 
—0.50 MeV -1 . Analyses that renormalize theoretically calculated 6 factors must include an additional systematic 
error in order to account for uncertainties in the true shape of the 6 factor when extrapolating to zero energy. It 
would be desirable to use the data themselves to determine the shape of the 6 factor. 

In this paper, we evaluate the modern prompt capture 7 ray and induced 7 Be activity measurements of the 
3 He(a,7) 7 Be cross section in a nearly model-independent way. We use well known physics to constrain the en- 
ergy dependence of the low energy cross section and disuss the treatment of systematic errors and their propagation 
into the final uncertainties, providing reliable values for 634 at energies relevant to solar neutrino production and 
BBN with realistic uncertainties derived from Markov Chain Monte Carlo calculations. 



II. MODERN DATA 



In the past 50 years, there have been many experimental efforts to measure the 3 He(a, 7) 7 Be cross section, including 
those of Holmgren and Johnston [32[ , Parker and Kavanagh [33l| . N agatani et al. [3j| , Krawinkel et al. [HI , Robertson 
et al. f36l . Volk et al. [13], Alexander et al. [H|, Osborne et al [39l|. and Hilgemeier et al. [i(| ■ With the exception of 
Ref. |32| . these measurements were all considered in the Adelberger et al. review Q, which recommended 6*34(0) = 
0.53 ± 0.05 keV b. This review noted a 2.5tr discrepancy between the weighted averages of the measurements based 
on prompt 7 ray detection and those based on induced 7 Be activity. At the time of this evaluation, we have the 
benefit of several subsequent, independent 7 Be activity measurements as well as two prompt 7 ray detection studies. 
These new measurements employed improved detectors, better background suppression, and underground accelerators. 
Moreover, they were published almost two decades after the last of the measurements considered in Ref. [f|. We opt 
therefore to restrict this analysis to modern experimental data, namely the recent measurements of Bemmerer et al. 
[HI, Brown et al. [H|, Confortola et al. [H|, Gyiirky et al. 0, and Nara Singh et al. 0. As the data of Ref. [U 
are included in the data of Ref. [3], we will refer henceforth only to the latter. The data considered here are shown 
in Table [J 



III. EVALUATION 



Several papers have discussed the energy dependence of nonresonant low energy radiative capture cross sections 
[2ll [3II l46j . When reactions are dominated by external capture, the astrophysical 6 factor exhibits a subthreshold 
pole of the form 

S(E) oc _L_, (1) 

where E is the relative energy and Q is the energy released by the reaction. This pole affects the convergence of the 
usual Maclaurin series expansions and the shape of cross sections at low energy, resulting in an upturn as E — > 0. If one 
wants better convergence, this subthreshold pole must be taken into account, expanding the quantity (E + Q)S(E). 

As shown in Ref. (47| . a completely model independent approach such as a Maclaurin series expansion does not 
work well when there are large discrepancies among data sets. In such cases one requires a physical constraint on 
the shape of the 6 factor. Since there have been many studies of the 3 He(a,7) 7 Be reaction, we can constrain the 
form of the fitting function beyond just a Maclaurin series using known physics. For instance, this radiative capture 
is dominated by El transitions at low energies [22j. Therefore, given the spin of nuclei involved, only I = and I = 2 
incoming partial waves contribute significantly to the external capture. 
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TABLE I: The modern data used in this evaluation. Shown are the S factors for capture into the ground (So) and first excited 
state (Si) individually, or the sum of the two contributions (Stot)- Relative systematic errors for prompt measurements (S p ), 
activity measurements (S a ), and those that are common to the two methods (S c ) are shown separately. 
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FIG. 1: Theoretical S34 calculations by Tombrello and Parker [22I ] (solid) and by Kajino [27I ] (dashed) are plotted relative to 
their values at E = in the upper graph. The deviations relative to the Kajino calculation are plotted in the lower graph. 
The dotted lines delineate the uncertainty due to the ± 6% theoretical uncertainty in the zero-energy S factor and logarithmic 
derivative estimated in Ref. [25| . 



Following Ref. [481 ] . we expand the remaining terms accounting for Coulomb and pole effects and find 

S ( E ) = 7T%^ Ml + aE + • • • ) 2 + a a (l + 4n 2 E/E G )(l + 16t: 2 E/E g )(1 + cE + ■ ■ -) 2 } . (2) 
E + Q 

Here, s an d S2 are the amplitudes of the s and d wave El capture components respectively. The terms with the 
coefficients a and c are higher order contributions. Keeping the high energy behavior of these terms equivalent (oc E 2 
to lowest order) demands c = if a is finite. Eq is the Gamow energy given numerically by Eq = 0.97913Z 2 Z 2 A 
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MeV where Z\e and Z^e are the charges of the reactants and A their reduced mass in atomic mass units. For S34, 
E G = 26.9437 MeV. 

The 3 He(a,7) 7 Be reaction can proceed via capture into either of the spin-orbit partners, the 7 Be ground state 
(Q = 1.5861 MeV) and the first excited state at 429 keV (Q = 1.1570 MeV). When determining best fits and zero- 
energy extrapolations, one must bear this in mind and use experiments sensitive to not only the total capture cross 
section, but also the partial capture cross sections into the ground and first excited states. We are therefore left with 
3 parameters for the S factor fits for each capture (i.e. ground or excited state). 

In general, modern experiments are dominated by systematic errors. As such, standard statistical treatments used 
to combine data from d ifferen t, possibly discrepant experiments are not valid. The traditional method of scaling the 
error in the mean by y^Jv, where v is the number of de gree s of freedom, has been shown to poorly estimate the 
error when systematics dominate the error budget [l5|, l49l l50j|. Such treatments lead to an underestimation of the 
true uncertainties. We require a more reliable prescription for this analysis. 

If one assumes that the dominant systematic error is the overall normalization error, several techniques can be used 
to estimate the true uncertainty. It was shown in Ref. flsT ] that a simple robust measure of systematic uncertainty is 
determined by the weighted dispersion of the fit. In another approach, d'Agostini [49| suggests floating the normaliza- 
tions of individual experiments weighted by their uncertainties as a way of compensating for the shortcomings of the 
standard statistical approach. This method works well, and agrees with the error estimate of Ref. [l5| . In the modern 
era, the size of systematic errors is often greater than or equal to the size of statistical errors. Hence systematic errors 
must be treated properly in order to obtain a reliable central value and uncertainties . 

We use a x 2 minimization procedure to determine the parameters of our calculated S34 that best fit the data. We 
break up the x 2 m t° two components, x 2 data and Xnorm' sucn that % 2 = x\ata + Xnorm- Thus for a single data set n, 

xl» = w ^-^ y and (3) 

Xnormi n ) = ( g ^ • (4) 

Here, a n and d n are the floating renormalization factor and the normalization uncertainty of the dataset n, S{Ei) is 
the theoretically calculated S factor at Ei , Si is the measured S factor at Ei , and Oi is the standard deviation of the 
measured S factor at Ei. 

In the case where two experiments have correlated normalizations, Xnorm ^ s modified and becomes 



2 , 1 9 s _ - l) 2 - 2,5 e 2 K 1 - 1)K 2 - 1) + S nl (an2 - l) 2 ] 

Xnorm \ n *-> r ^ 2 ^2 _ 



where S n i and 5 n 2 are the total systematic errors of the data sets nl and nl and S c is the systematic error common 
to both data sets. This equation is the result of inverting the covariance matrix. 

We use a Markov Chain Monte Carlo (MCMC) approach to evaluate the data, following the two-step procedure 
laid out by Ando et al. [5lJ. First, we find the best fit by varying the parameters randomly, with ever decreasing step 
sizes until the results converge to a predetermined number of significant figures. Then, with step sizes determine d by 
Ax 2 = 1 variations, we start a random walk away from the best fit point using the Metropolis algorithm [52L l53|. 
Once the sample variance of each of the parameters converges within the specified resolution, we stop the MCMC. 
We tested convergence with chains of length up to 10 8 , finding convergence after about 10 6 steps. We use the 10 8 
step chain for our final results. 



IV. RESULTS 



To test the validity of the adopted three parameter S 34 model, we found the best fits for both two and four 
parameter models as well, fixing a — for the two parameter fit and letting both a and c vary for the four parameter 
fit. We find best fits with Xtot = 85,63,61 or x 2 otl v — 2.74,2.17,2.26 for the two, three, and four parameter fits, 
respectively. One can see readily that given the modern data considered here, we can not adequately determine more 
than 3 parameters. Also, the two parameter fit is clearly unrepresentative of the data, whilst the quality of the four 
parameter fit is no better than that of the three parameter fit. Hence including a fourth parameter is unwarranted. 
The uniqueness of the solution was tested by starting the minimization procedure from different parameter values, 
finding the same minimum. 
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Parameter 


Mode Mean ± Std. Dev. 


Norm. Error 


a p (Brown) 
a a (Brown) 
q p (Confortola) 
a a (Confortola) 
Q a (Gyiirky) 
a a (Nara Singh) 


0.95 0.95 ±0.02 
0.95 0.95 ±0.02 
0.99 0.99 ±0.02 

1.01 1.01 ±0.02 

1.02 1.02 ±0.02 
1.04 1.04 ±0.02 


0.0387 

0.030 

0.038 

0.032 

0.031 

0.0363 


so(gs) 

82 (gs) 

o(gs) 


0.406 0.406 ±0.009 
0.007 0.007 ±0.001 
-0.207 -0.203 ±0.038 




s (ex) 
s 2 (ex) 
o(ex) 


0.163 0.163 ±0.004 
0.004 0.004 ±0.001 
-0.134 -0.127 ±0.055 




Xnorm 
2 

Xdata 


4.83 6.48 ±2.02 
58.38 69.06 ±4.70 


Nnorm — 6 
N data = 41 


Xtot 


63.21 75.54 ±4.95 


Npar = 6 



TABLE II: Results of the Markov Chain Monte Carlo parameter estimation. 



After the best fit is found by minimizing \ 2 an d the MCMC has converged, we find the limits of central 
confidence intervals by determining parameter values that obey the relation P(v/2, Ax 2 /2) = CL (54|. Here 
P(a,x) = 7(a, x)/T(a) = 1 — T(a, x)/T(a) is the regularized incomplete gamma function [55|, v = N par + N norm 
is the total number of varied parameters and normalizations, A% 2 = "x 2 — Xmiw ano - CL is the desired confidence 
level. 

The results of this analysis are presented in Table II and Figure [2] We find a mode and central 68.3% CL 
interval of S" 34 (0) = 0.580 ± 0.043 keV b (AS 34 (0) / S 34 (0) = 7.4%). We find a mean and standard deviation of 
534(0) = 0.580 ± 0.013 keV b (AS , 34 (0)/ '5 34 (0) = 2.2%). The small size of the interval given by the mean ± the 
standard deviation compared to the central 68.3% CL results from marginalizing over all parameter distributions. 
The standard deviation here corresponds to defining errors with A% 2 = 1, known to underestimate uncertainties when 
there is more than 1 degree of freedom. The central 95.4% confidence interval is 6*34(0) = 0.580 ± 0.054 keV b. The 
size of this relative uncertainty is slightly smaller than that given in the evaluation of Adelberger et al. [6| , which was 
a la error. Hence the modern data considered here permit a considerably more precise recommendation for £34 (0), 
even when allowing the shape of the S factor to be determined by the data rather than a theoretical model and when 
fitting the branching ratio as well as the total S factor. 

Using the best fit parameters and the definitions in [la ], we find a discrepancy error of 5.3% and a normalization 
error of 3.0%, yielding a total systematic uncertainty of 6.1%, in very good agreement with the present MCMC results, 
considering that this error does not take into account discrepancies in the branching ratio Si/ So, as the MCMC method 
does. The MCMC results for the branching ratio are shown in Figure (3) The s- and d-wave contributions to the 
ground and excited state transitions determined in the MCMC analysis are shown in Figured) 

We have calculated the thermally averaged rate of the 3 He(a,7) 7 Be reaction per particle pair A34 and fit it using 
the form 



, . ^2/3 . ^,4/3 

1 + mTg' +n 2 T g ' 



a ± a 2 /T 9 1/J ± a 6 lnT 9 -) — Tyr^cm 3 moL 1 b" 1 . (6) 

1 ± drTf 3 + d 2 T 9 4/3 



A 34 = exp 

The fit parameters are shown in Table III. This functional form provides a better fit to the reaction rate than the 
fitting functions of Caughlan and Fowler [5(| and REACLIB [U HH and agrees with the numerically calculated rate 
within 0.5%. 

Figure [5] compares the energy dependence of S34 determined here from a nearly model-independent analysis of 



modern data with that calculated by Kajino et al. 27J. It is apparent that in the energy range of principal interest 
below 500 keV, the shapes of the two curves differ substantially. Extrapolation using the RGM model of Kajino from 
the energy range important for BBN to that of interest in the Sun would differ by some 5% from extrapolation based 
on the S factor determined here. This is quite consistent with the estimates of the theoretical error given in Ref. [25| . 
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Parameter 


Low 


Adopted 


High 


fflo 


15.531721 


15.609867 


15.679639 


m 


-0.100208 


-0.020478 


0.037757 


n 2 


0.235187 


0.211995 


0.196645 


di 


0.114322 


0.255059 


0.353050 


d 2 


0.373802 


0.338573 


0.316019 



TABLE III: Shown are the fit parameters for the recommended thermal rate (A34) and its central 68.3% confidence level limits, 
with a 2 = -12.82707707, a e = -2/3, and a = In [1.03762 x lO^WkcV^bT 1 ] fixed. 
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FIG. 2: (Color online) The best fit S factor for the 3 He(a,7) 7 Be reaction (solid curve). Also shown are the central 68.3% 
confidence level limits (dashed curves). The experimental data and their total uncertainties are plotted (Activity-open points, 
Prompt-solid points). 
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FIG. 3: (Color online) The best fit branching ratio Si/ So for the 3 He(a,7) 7 Be reaction (solid curve). Also shown are the central 
68.3% confidence level limits (dashed curves). The experimental data and their total uncertainties are plotted. 



10 




0.5 1 1.5 
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FIG. 4: (Color online) This figure shows the fitted contributions of various partial waves to the S factor. The s-wave (short 
dashed), d-wave (dotted) and total (solid) S factors for capture into the ground state and the s-wave (long dash-dotted), d-wave 
(short dash-dotted) and total (long dashed) S factors for capture into the first excited state are all shown. 

V. CONCLUSIONS 

We have performed a minimally model dependent analysis of modern 3 He(a,7) 7 Be data. This analysis properly takes 
into account the systematic errors of potentially discrepant data and yields reliable central values and uncertainties. 
At the 68.3% confidence level, we find 5 34 (0) = 0.580 ±0.043 keV b and S'(0)/S(0) = -0.92 ±0.18 MeV" 1 . We have 
used this value of 634(0) to compute the thermally averaged rate of the 3 He(ev,7) 7 Be reaction per particle pair, and 
have fit this rate with an analytical form accurate within 0.5%. 

This reaction rate is about 9% higher than the recommendation of [6j, but is quite consistent with it. It is also 
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FIG. 5: (Color online) Comparison of the shape of the RGM calculation of Kajino [27j ] with the present structure model- 
independent analysis. 

perfectly consistent with the recommendations of References [H^, |60j]. Moreover, access to precise modern data has 
allowed us to minimize the uncertainties associated with extrapolation and arrive at a more precise and better founded 
value that does not depend on nuclear structure models. 

Our recommended reaction rate implies that the predicted 7 Be and 8 B solar neutrino fluxes should both be increased 
by 8% compared with predictions based on the 534(0) value from Ref. Q. Similarly, the present analysis implies that 
the value of 534 in the energy range relevant to big bang nucleosynthesis is 17% larger than the value adopted in Ref. 
(TH . resulting in a 17% increase in the predicted primordial 7 Li abundance. The uncertainties found here reduce the 
total error in the 7 Li abundance prediction by a factor of ~ 2. The discrepancy between the predicted and observed 
primordial 7 Li abundances increases from ~ 3<r to ~ 5a [6l|, implying a rather serious disagreement that must be 
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resolved. 

In order to improve our knowledge of the rate of the 3 Hc(a,7) 7 Be reaction, future experiments should collect data 
at at least six different energies. No modern data exist in the Gamow window for big bang nucleosynthesis, so this 
would be of primary interest. Information on capture into the ground and first excited states would also be useful, 
arguing in favour of prompt 7 ray measurements rather than those based on induced 7 Be activity. Finally, precise 
measurements from 1-3 MeV would allow better constraints on the shape of the S factor at the low energies of 
astrophysical interest. 
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